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SUMMARY 


A  procedure  is  presented  for  the  efficient  computational  solution  of 
linear  programs  having  a  certain  structured  property  characteristic  of  a 
large  class  of  problems  of  practical  interest.  This  property  makes  possible 
the  decomposition  of  the  problem  into  a  sequence  of  small  linear  programs 
whose  Iterated  solutions  solve  the  given  problem  through  a  generalization 
of  the  simplex  method  for  linear  programming. 
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THE  DECOMPOSITION  ALGORITHM  FOR  LINEAR  PROGRAMMING^ 

1.  THE  DECOMPOSED  LINEAR  PROGRAM 

Many  linear  programming  problems  of  practical  interest  have  the 
property  of  being  composed  of  separate  linear  programming  problems,  which 
are  tied  together  by  a  number  of  constraints  that  is  considerably  less  than 
the  total  number  imposed  on  the  problem.  When  the  suitably-ordered  matrix 
of  coefficients  of  such  a  problem  is  displayed  in  the  usual  way,  a  pattern 
emerges  like  that  shown  in  Fig.  1.  In  this  figure  the  constraint  matrix  has 


B 

n 


b 

n 


Fig.  1  —  The  original  problem 


^A  preliminary  version  of  this  research  memorandum  appeared 
originally  under  the  title  A  Decomposition  Principle  for  Linear  Programs, 
The  RAND  Corporation,  Paper  P-1544,  November  10,  1958,  and  was 
presented  at  The  RAND  Symposium  on  Mathematical  Programming  in 
March,  1959. 
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been  partitioned  into  non— zero  blocks  A  and  B  ;  the  "right-hand  side" 

J  J 

column  of  constants,  b,  b^,  b^,  is  arrayed  correspondingly;  and 

the  "costs"  ■ —  the  coefficients  of  the  objective  form  —  are  partitioned 
into  the  blocks  c^,  c^, 


• » 


n 


More  particularly,  suppose  that  we  have  the  following; 


an  m-vector  b. 


A. ,  an  m  by  n.  matrix, 

B. ,  an  m.byn.  matrix, 

]  3  J  ’ 

c.,  an  n,— vector,  and 
3  3 

X.,  a  variable  n.— vector, 
3  3 


for  each  j  =  1,  .  .  . ,  n. 


Then  the  problem  illustrated  in  Fig.  1  is  posed  as  a  linear  programming 

problem  in  E.n,  variables,  subject  to  m  +  E,m.  constraints.  It  can  be 
3  3  3  3 

summarized  as  follows: 

The  Decomposed  Program.  Find  the  vectors  x.  (j  =  1,  .  .  . ,  n) 

3 

satisfying  the  constraints 


(1)  E.A.x.  =  b  and  x.  >  0  (for  all  j), 

3  3  3  3  - 

(2)  B.x,  =  b,  (for  all  j), 

3  3  3 

which  minimize  the  linear  form 


(3) 


E . c ,x. . 
j  3  3 
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It  would  seem  that  each  of  the  n  sets  of  constraints  of  (2)  constitutes 
a  Bubproblem  of  secondary  importance  to  the  whole  program,  and  that  these 
sets  should  be  studied  mainly  through  the  restrictions  they  impose  on  the 
activities  of  the  joint  constraints  (1).  Pursuing  this  point  of  view  leads  us 
from  the  decomposed  program  to  the  equivalent  extremal  problem  of  the 
next  section.  For  the  sake  of  simplicity  in  formulating  the  extremal  problem, 
it  is  assumed  that 

<4)  Sj  .  I  X.  >  0,  B.X.  =  bj  ) 

is  bounded  for  each  j.  Only  minor  changes,  which  will  be  discussed  in 
Sec.  4,  are  required  in  order  to  remove  this  restriction  in  the  treatment 
which  follows. 

The  results  of  this  paper  arise  from  the  statement  of  the  decomposed 

program  as  a  "generalized  linear  programming  problem"  in  n  variables, 

in  which  the  m-element  column  of  coefficients  associated  with  each  variable 

is  drawn  freely  from  the  convex  set  S,  instead  of  being  fixed  as  in  the 

J 

ordinary  linear  programming  problem.  This  generalized  problem  is 
developed  more  fully  elsewhere  [5]  . 

2.  THE  EXTREMAL  PROBLEM 

In  this  section  we  formulate  a  type  of  linear  programming  problem 
in  which  those  features  of  the  decomposed  program  (1-3)  that  pertain  to  its 
joint  constraints  are  brought  to  the  fore.  In  the  new  problem,  the  extreme 
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points  of  the  sets  defined  by  the  relations  (2)  will  yield  the  data  for  the 
problem. 

For  j  ■  1,  .  . . ,  let  . ""jK  ^ 

«} 

extreme  points  of  the  convex  polyhedron  defined  by  the  conditions 

X,  >  0,  Fi.x,  “  b.  ;  also  define 
j  -  J  J  3 


(5) 


P..  =  A.x 
]k  j  3k 

c.,  =  c.x 
]k  3  3k 


V  for  k  =  1,  .  .  . ,  K 


2.  1.  The  Extremal  Program 

Find  numbers  s (j  °  1,  •  •  •  >  ^  *  *  ’  >  satis^'ving 


(6) 

and 

(7) 


L  P  s.  =  b,  s,,  >  0  (for  all  j,  k), 

3,k  jk  3k  jk  - 


(for  all  j), 


which  minimize  the  linear  form 


(8)  E.  ,  c.,  s.,  . 

3,k  3k  3k 

The  relation  of  the  extremal  problem  to  the  original  problem  lies  in 

the  fact  that,  since  S.  is  a  bounded,  convex  polyhedral  set,  any  point  x^ 

3 

of  S  may  be  written  as  a  convex  combination  of  its  extreme  points,  that 
3 

is,  as  where  {s.^,  satisfy  (7);  and  the  expressions 

(6)  and  (8)  are  just  the  expressions  (1)  and  (3)  of  the  decomposed  problem 


Pri 
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Fig.  2  —  The  extremal  problem 
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constraints  of  the  original  problem  have  gone  over  into  the  m  constraints 

(6)  constituting  the  upper  block  in  Fig.  2,  and  the  constraints  of  the 
th 

j  subproblem  have  gone  over  into  single  constraints  of  the  form  (7).  The 
reduction  in  the  total  number  of  constraints  is  sizable  in  case  the  m.  are 

J 

large,  and  it  is  on  this  fact  that  the  computational  efficiency  of  the 

decomposition  algorithm  relies.  The  reduction  appears  to  have  been 

accomplished,  however,  by  greatly  enlarging  the  number  of  variables  in 

the  problem  from  the  original  E.n.  to  the  number  E  K.  ;  the  proposed 

J  J  J  J 

method  would  be  of  little  Interest  or  value  if  it  were  not  possible  effectively 
to  reduce  this  number. 


3.  THE  DECOMPOSITION  ALGORITHM 

The  central  idea  of  the  decomposition  algorithm  is  that  the  extremal 

linear  programming  problem  of  (6—8)  and  Fig.  2  can  be  solved  by  the 

simplex  method  for  linear  programming  without  prior  calculation  of  all 

the  data  given  in  the  statement  of  the  problem.  Typically,  only  a  small 

number  of  the  E.K.  variables  present  will  ever  play  an  active  role  in  the 
J  J 

course  of  solving  the  problem,  and  as  will  be  seen,  the  coefficients  needed 
in  handling  just  these  variables  may  be  generated  at  the  time  they  are  to 
come  under  consideration.  The  principle  of  coefficient  generation  was 
first  employed  in  the  study  of  multi— commodity  network  flows  [3]  . 

Let  us  repeat  here  the  essential  features  of  the  simplex  method  that 
will  be  used.  The  phenomenon  of  degeneracy  plays  the  same  role  in  the 
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extremal  problem  that  it  does  in  any  linear  programming  problem,  and 
standard  methods  [l]  can  be  invoked  to  handle  it  when  necessary,  so  that 
it  need  not  be  discussed  here.  The  problem  discussed  in  the  next  section  — 
that  of  finding  a  first  feasible  solution  to  the  extremal  problem  —  is  similar 
to  the  same  problem  for  a  general  linear  program.  It  will  be  assumed 
that  the  usual  data  are  at  hand  for  performing  one  typical  step  of  the  simplex 
method  on  the  extremal  problem. 

Since  the  extremal  problem  has  m  +  n  equation  constraints,  there 
will  be  at  hand  m  +  n  columns  { (m  +  n)— component  vectors  corresponding 
to  extreme  points)  that  constitute  a  "feasible  basis."  These  columns  are 
linearly  independent,  and  the  unique  solution  of  the  constraint  equations 
(6,  7),  which  is  obtained  by  setting  to  zero  those  variables  associated  with 
all  other  columns,  is  nonnegative.  If  the  revised  simplex  method  is  used 
in  performing  the  calculations,  there  will  also  be  at  hand  the  (m  +  n)-vector 
of  "prices"  (tt;  n)  — the  m-^ector  tt  being  associated  with  the  first  m 
constraints  and  the  n-vector  tt  with  the  remaining  n,  as  in  Fig.  2.  In  the 
general  linear  program,  the  inner  product  of  the  price  vector  with  any 

column  of  the  basis  must  equal  the  cost  associated  with  that  column.  In 

““  th  " 

the  case  of  the  extremal  problem,  if  we  let  be  the  j  component  of  ir , 

this  relationship  can  be  written  as 


(10) 


TT  P.,  +  TT,  =  C  .. 

Jk  3  3K 


th 

for  basic  columns  drawn  from  the  j  partition,  3=1,  •  •  • ,  n 
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One  step  of  the  simplex  method  iteration  for  solving  the  extremal 
problem  would  be  performed  as  follows:  Find  a  column  of  the  constraint 
matrix  whose  "reduced  cost"  is  negative,  that  is,  for  which 


(11) 


0  . 


(Commonly,  the  column  which  minimizes  the  reduced  cost  is  chosen. ) 
Add  this  column  to  the  current  basis,  and  delete  one  column  from  the 
basis  in  such  a  way  that  the  new  basis  is  still  feasible.  If  no  column 
gatisfying  (11)  can  be  found,  then  the  current  solution  solves  the 

extremal  problem.  Otherwise,  the  simplex  method  gives  the  appropriate 
rules  for  the  removal  of  a  column  from  the  basis,  and  for  the  calculation 
of  the  new  prices  (tt;  ^)  associated  with  the  new  basis  with  which  the  next 
iterative  step  can  begin. 

A  procedure  for  applying  the  simplex  algorithm  to  the  extremal 
problem  without  having  all  the  data  of  that  problem  at  hand  can  now  be 
stated.  It  is  supposed  that  only  the  m  +  n  columns  of  a  feasible  basis 
for  the  problem  are  given. 


3.  1.  The  Decomposition  Algorithm:  Iterative  Step 

Given  m  +  n  columns  of  the  form  •  •  •  >  •  •  •  >  0) 

constitute  a  feasible  basis  for  the  extremal  problem,  their  associated 


costs  c.,  ,  and  prices  (tt;  tt)  satisfying  (10);  for  each  j  -  1,  .  .  .,  n,  let 
3k» 

X.  be  an  extreme  point  of  S.  minimizing  the  linear  form 
3  J 
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(12) 


(c.  —  7rAj)x^,  under  the  constraints  >  0,  B^x^  *  ; 


and  let  x  be  such  that 


(13) 


(c  -ttA.  )x.  -TT  =  Min  (c  -  77- A,)x  -  TT  . 
Jq  ^0  ^0  ^0  J  L  j  J  J  JJ 


If  6  <  0,  form  the  new  column  and  its  associated  cost  for  the  extremal 
problem  as 


(14)  (A,  X  ;  0,  •  ‘  ‘  f  Ij  •  •  •  »  0)  > 

Jq  h  •’o  Jq 

add  this  column  to  the  basis,  and  form  a  new  basis  and  new  prices  using 
the  rules  of  the  simplex  method.  If  6  >  0,  terminate  the  algorithm; 

{s  }  solves  the  extremal  problem,  and  the  relations  (9)  give  the  solution 

jk 

of  the  original  problem. 


Theorem.  The  decomposition  algorithm  terminates  in  a  finite 
number  of  iterations,  yielding  a  solution  of  the  extremal  problem. 

Proof.  Since  the  termination  of  the  simplex  method  as  applied  to 
the  general  linear  programming  problem  in  a  finite  number  of  iterations 
is  known  [ij  ,  it  is  sufficient  to  show  that  the  rules  of  the  iterative  step 
of  the  decomposition  algorithm  yield  a  column  satisfying  the  criterion  (11), 
if  this  is  possible.  For  the  column  defined  by  (14),  the  left-hand  side  of 
(11)  is  simply  c.  x.  -ttA,  -  tt  .  =6,  which  is  as  small  as  possible. 

^The  lemma  of  Sec.  2  shows  that  a  solution  of  the  original  problem  is 
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obtained  when  a  solution  of  the  extremal  problem  has  been  found  by  this 
algorithm. ) 

In  summary,  a  cycle  of  the  decomposition  algorithm  can  be  stated 
this  way: 

Given  m  +  n  columns  constituting  a  feasible  basis  for  the  extremal 
problem,  use  the  prices  associated  with  that  basis  to  form  the  modified 


costs,  c^  -  TT  Aj,  for  each  of  the  subproblems. 

J  J 

Minimize  (c  -7rA.)x.  for  x  >  0,  B  x  »  b  .  If  6  <  0,  so  that 
J  3  3  3  3  3  3 

the  original  problem  is  not  solved,  then  from  an  appropriate  subproblem 


solution,  construct  a  new  column  for  the  extremal  problem  and  form  a 
new  feasible  basis  incorporating  that  column,  deleting  another  column 
from  the  basis  in  accordance  with  the  rules  of  the  simplex  method. 


4.  DETAILS  ON  USE  OF  THE  ALGORITHM 

This  section  considers  some  details  regarding  use  of  the  decomposition 
algorithm:  how  to  get  the  algorithm  started;  how  to  deal  with  unbounded 
solutions  of  the  subproblems,  which  may  occur  when  the  boundedness 
restriction  is  removed;  some  variations  of  the  selection  technique  (12) 
and  of  methods  for  decomposing  a  problem. 

4.1.  Initiating  the  Algorithm 

The  decomposition  algorithm  can  be  started  with  precisely  the  same 
device,  called  Phase  One,  that  is  used  for  the  general  linear  programming 
problem.  This  device  consists  in  augmenting  the  problem  with  m  +  n 
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"artificial"  variables  in  terms  of  which  an  initial  basic  feasible  solution, 
and  the  prices  associated  with  the  corresponding  initial  feasible  basis,  are 
readily  found.  The  decomposition  algorithm  can  then  be  applied  to  the 
problem  of  removing  the  artificial  variables.  After  this  has  been  done, 
the  required  Initial  conditions  for  the  ordinary  application  of  the  algorithm 
are  automatically  met. 

For  i  =  1,  . .  . ,  m  +  n,  let  y^  be  a  nonnegative  variable;  let  1^  be 
the  l^'^  column  of  the  (m  +  n)— order  identity  matrix;  and  let  c^^  =  1  be  the 
cost  associated  with  the  variable  y^  .  For  Phase  One,  replace  all  the  costs 
{c  }  of  the  original  problem  with  zero  vectors.  It  can  be  assumed  without 

J 

loss  of  generality  that  the  right-hand  side  vectors  b,  {b^}  are  all  non¬ 
negative. 

Designating  {I^}  as  the  initial  feasible  basis,  employ  the  decomposition 
algorithm  in  the  minimization  of  the  linear  form  [Note  that  the  Initial 

feasible  solution  and  the  initial  prices  are  given  by  (y^,  .  . . ,  “ 

(b;  bj,  .  .  . ,  b^)  and  {tt;  tt)  =  (1,  .  .  . ,  1).  ] 

If  the  form  ^^y^  cannot  be  reduced  to  zero,  then  the  extremal  problem 
has  no  feasible  solution;  if  the  form  can  be  reduced  to  zero,  then  a  feasible 
solution  is  a  fortiori  at  hand.  Typically,  this  Phase  One  will  have  ended 
with  none  of  the  artificial  columns  left  in  the  basis.  (It  occasionally  happens 
that  an  artificial  column  remains  in  the  basis  although  its  variable  is,  of 


course,  zero.  Handling  this  case  requires  the  use  of  an  additional  constraint, 
described  elsewhere  [l]  ,  in  the  problem.)  The  cost  vectors  of  the  original 
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problem  are  then  restored,  and  Phase  Two,  the  application  of  the  algorithm 
to  the  proper  extremal  problem,  can  proceed  as  in  Sec.  3. 

4.  2.  Extension  to  the  Unbounded  Case 

It  was  assumed  in  (4)  that  the  constraint  set  of  each  subproblem 
was  bounded.  If  this  is  not  the  case,  then  it  may  happen  in  applying  the 
decomposition  algorithm  that  an  unbounded  solution  is  obtained  for  the 
subproblem: 

(15)  Minimize  (c^-ttA^Ix^  under  the  constraints  x^  >  0,  B^x^  =•  b^  . 

The  only  change  incurred  in  the  algorithm  by  this  is  that  the  rule  for  forming 
a  new  column  for  the  extremal  problem  must  be  extended.  An  unbounded 
solution  of  the  problem  (15)  will  cause  the  simplex  method  to  yield  a  vector 
V.  such  that 

(16)  y.  >  0,  B.y^  =  0,  and  <  0- 

The  relationships  (16)  give  the  direction  of  an  infinite  ray  of  feasible 
solutions  along  which  the  cost  form  of  (15)  proceeds  to  -  oo  . 

In  this  case,  the  column  and  cost  to  be  added  to  the  extremal  problem 

have  the  form 

(17)  (A  y  ;  0,  .  . . ,  0)  and  c  y . , 

J  V 

instead  of  the  form  (14).  The  "l”  has  been  omitted  from  the  column  so 
that  the  variable  s.,  associated  with  this  column  will  not  be  constrained 
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bv  the  relation  L,  =■  1,  imposed  upon  the  columns  created  from  bounded 
solutions.  Any  nonnegative  multiple  of  this  column  is  admissible,  corre¬ 
sponding  to  the  fact  that  any  nonnegative  multiple  of  the  vector  y^  is 
admissible  in  the  representation  of  the  constraint  set  of  (15)  as  the  sum 
of  a  bounded  polyhedron  and  a  polyhedral  convex  cone.  It  is  easy  to  check 
that  with  this  extra  freedom  the  decomposition  algorithm  still  consists 
only  of  applying  the  simplex  method  to  the  extremal  problem  thus  extended, 
and  hence  the  algorithm  has  the  same  termination  properties  as  in  the 
bounded  case,  since  the  number  of  possible  rays  generated  by  the  simplex 
solutions  of  the  subproblems  is  finite. 


4.  3.  Variations  of  the  Decomposition  Process 

It  is  possible  to  pursue  the  decomposition  algorithm  in  different  ways 
for  the  sake  of  computational  efficiency.  One  suggestion  is  easily  made: 


retain  as  many  of  the  solutions  x^  to  the  subproblems  (12)  as  would  provide 
new  basic  columns  for  the  extremal  problem;  that  is,  for  each  j  such  that 


form  the  column  (A.x;  0,  .  .  . ,  1,  •  •  •  ,  0)  and  its  cost  c  x  ,  and  adjoin 

.1  j 

all  these  to  the  given  basic  columns,  employing  the  simplex  method  to  find 
a  minimum  among  this  extended  set  of  columns.  This  procedure,  while 
requiring  perhaps  more  simplex-method  iterations,  should  require  fewer 
subprogram  solutions  to  be  found. 
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It  should  be  noted  that  finding  the  n  subproblem  solutions  of  (12)  at 
each  iteration  of  this  procedure  may  not  require  the  complete  solution  of 
n  linear  programs.  If  the  prices  in  the  extremal  problem  do  not  change 
much  from  one  iteration  to  the  next,  then  the  new  subproblem  solutions 
may  not  differ  much,  if  at  all,  from  their  former  ones;  in  that  case,  the 
new  solutions  can  be  found  with  little  labor.  Also,  in  many  practical  cases, 
the  subproblems  themselves  have  such  structure  that  their  solutions  are 
readily  found  by  means  of  special  devices  —  for  example,  when  the  sub- 
problems  are  transportation  problems. 

Another  variation  in  the  use  of  decomposition  depends  on  the  fact 
that  a  given  linear  programming  problem  may  be  decomposed  to  various 
degrees.  For  example,  in  the  original  problem  (1-3)  pictured  in  Fig.  1, 
the  first  two  partitions  might  have  been  considered  collectively  as  one, 
embracing  the  single  vector  (x^^;  x^).  Little  would  be  changed  in  the 
description  of  the  decomposition  algorithm,  except  that  the  first  subproblem 
of  the  family  (12)  would  read  as  follows: 

(18)  Minimize  (c ^  -  tt  Aj)Xj  +  (c2  A2)x2  under  the 

restrictions  (x^^,  x^)  >  0,  ~  ^2^2  "  ^2' 

the  column  constructed  therefrom  would  be 

(19)  (A^x^  +  A^x^;  0,  ...,  1,  ...,  0),  where  (x^,  x^)  solves  (18). 
The  two  variables  of  the  problem  (18),  however,  are  completely  independent. 
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so  that  its  solution  (x^;  x^)  is  psi  composed  of  the  solutions  x^,  x^,  of 
the  two  separate  problems  of  the  ordinary  form  for  j  “  1,  2.  The  result, 
then,  of  "aggregating”  the  first  two  partitions  has  been  to  aggregate  their 

resultant  columns  according  to  the  formula  (19). 

If  this  aggregation  is  carried  as  far  as  possible,  all  the  n  subproblems 

may  be  taken  together  as  one.  On  account  of  the  independence  of  the  n  parts 
of  this  one  subproblem,  finding  the  subproblem  solution  involves  the  same 
work  as  before.  An  advantage,  however,  lies  in  the  fact  that  the  extremal 
problem  now  has  only  m  +  1  constraints,  instead  of  m  +  n,  and  thus  is 
easier  to  handle.  The  decomposition  of  a  linear  programming  problem 
into  a  problem  having  only  one  partition,  when  applied  to  the  classical 
transportation  problem,  yields  an  algorithm  which  promises  considerable 
efficiency  in  the  case  of  transportation  from  a  small  number  of  origins  to 


large  number  of  destinations  [4]  . 

Finally,  It  should  be  mentioned  that  the  linearity  of  the  cost  functions 


c.x  for  the  partitions  of  the  problem  is  not  essential  for  the  application 
of  the  decomposition  algorithm.  Suppose  instead  that  in  each  partition  this 


function  were  replaced  by  the  function  Ux.).  which  need  only  be  convex 
for  success  of  the  method.  The  decomposition  algorithm  would  be  unchanged, 


except  that  the  subproblems  (12)  would  assume  the  form 


(20)  Minimise  t.(x.)  -  x  A.x,  under  the  constraints  x,  >  0,  B.x^  »  b^ . 
Solutions  X,  of  these  problems  would  then  be  used  to  construct  new  extremal 
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columns  by  formula  (14)  with  associated  costs  f,(x.).  In  the  nonlinear 
case,  termination  of  the  decomposition  algorithm  is  not  assured,  but 
convergence  to  a  solution  of  the  problem  is.  The  arguments  used  differ 
from  those  above;  this  extension  of  the  decomposition  algorithm  is  given 
elsewhere  [  5  ]  . 


5.  A  NUMERICAL  EXAMPLE 

The  small  problem  we  shall  solve  by  means  of  the  decomposltlpn 
algorithm  is  somewhat  artificial,  but  it  will  suffice  to  show  the  utility  of 
the  method  even  in  the  case  of  a  single  partition. 

A  homogeneous  commodity  is  to  be  transported  from  each  of  two 
sources  to  each  of  four  destinations.  The  costs  of  transporting  one  unit 
of  the  commodity  from  a  given  source  to  a  given  destination  appear  as 
entries  in  the  matrix  of  Fig.  3.  Beside  the  matrix  and  below  it  are  written. 


Fig.  3  —  Data  for  the  example 
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respectively,  the  amounts  of  the  commodity  available  from  each  source 
and  demanded  at  each  destination.  One  requirement  is  added  which  prevents 
this  from  being  an  ordinary  transportation  problem:  it  is  required  that  the 
sum  of  three  times  the  quantity  shipped  from  source  1  to  destination  3  and 
twice  the  quantity  shipped  from  source  2  to  destination  2  be  10. 

In  order  to  formulate  this  as  a  linear  programming  problem,  denote 


by  t^j  (i  *  1|  2;  j  =  1, 


4)  the  quantity  to  be  shipped  from  source  1  to 


destination  j .  The  coefficients  of  all  the  constraints  are  shown  in  Fig.  4 


Fig.  4  —  Decomposition  of  the  example 

since  there  le  only  one  partition,  subscripts  are  not  needed  to  describe  the 
parts  ot  the  matrix.  The  row  constituting  the  matrix  A  represents  the  Joint 
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constraint  on  t.^  and  t^g  •  The  block  B  consists  of  six  equations  asserting 
that  the  net  flow  from  each  source  must  be  just  the  amount  available,  and 
that  the  net  flow  to  each  destination  must  be  what  is  demanded. 

A  has  only  one  row,  and  the  problem  has  only  one  partition,  so  that 
m  ■  1  and  n  =*  1.  The  extremal  problem  is  thus  a  two-equation  problem, 
and  the  price  vector  (tt;  tt)  has  two  components.  Of  course,  “  x. 

The  single  subproblem  of  (12),  the  problem  of  minimizing  (c  —  n  A)x  under 
the  constraints  Bx  =  b,  becomes  here  just  that  of  solving  the  transportation 
problem  shown  in  Fig.  5,  since  the  subproblem  constraints  are  precisely 
those  of  the  transportation  problem  of  Fig.  3.  Rapid  special  methods  for 
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Fig.  5  —  The  subproblem  as  a  transportation  problem 


solving  transportation  problems  are  well  known  [2]  ;  we  need  not  go  into 
them  here. 

Beginning  with  a  basic  feasible  solution  for  the  extremal  problem 
constructed  from  the  artificial  columns  Ij,  (see  Sec.  4),  the  major 
data  generated  in  the  decomposition  algorithm  solution  of  this  problem 
are  given  in  the  table  below.  We  have  omitted  only  the  inverse  of  each 
feasible  basis,  which  is  generally  carried  along  in  simplex  method  calculations. 
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